Magnetic domain formation in itinerant metamagnets 
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We examine the effects of long-range dipolar forces on metamagnetic transitions and generalize 
the theory of Condon domains to the case of an itinerant electron system undergoing a first-order 
metamagnetic transition. We demonstrate that within a finite range of the applied field, dipolar 
interactions induce a spatial modulation of the magnetization in the form of stripes or bubbles. Our 
findings are consistent with recent observations in the bilayer ruthenate Sr3Ru207. 



Itinerant electron systems exhibiting metamagnetism, 
quantum criticality and non-Fermi-liquid behavior con- 
tinue to attract widespread experimental 0, 13- and 
theoretical 0, 0, B| interest with the bilayer ruthenate 
Sr3Ru2C>7 as the most prominent example. Particularly 
surprising was the recent discovery of a new phase with 
increased resistivity as the quantum critical point is ap- 
proached |3-Lj|. Different theories have been proposed for 
the nature of this new phase In this Letter, we show 
that dipolar magnetostatic forces in an itinerant meta- 
magnet lead to domain formation in a small yet finite 
region of the temperature field phase diagram. These 
domains differ in the magnitude rather than the direc- 
tion of the magnetization. Magnetic domain formation 
of this type was first proposed by Condon to explain de 
Haas - van Alphen measurements in Be []|. Much later, 
Condon domains, which are formed and destroyed in each 
de Haas - van Alphen cycle, were observed directly using 
muon-spin-rotation (/iSR) spectroscopy [2j and recently 
the domain pattern was observed in Ag by micro Hall 
probes |9j- The theory has been further refined in the 
context of de Haas - van Alphen effect [l0(, but to our 
knowledge, magnetic domain formation has not yet been 
considered in the context of metamagnetism. 

We first consider briefly the standard phenomenologi- 
cal theory of itinerant metamagnetism starting from the 
free energy per unit volume /(M, T) as a function of 
magnetization (produced by spin polarization of a metal- 
lic electron band) and temperature. The thermodynamic 
equation of state is fioH ln = 8m f, where fiQ is the vac- 
uum permeability and H ln is the magnetic field inside the 
sample. The differential susceptibility is x — dM/dH™, 
or equivalently HoX 1 — 9f jf. By definition, metamag- 



netism occurs if the susceptibility has a maximum at a fi- 
nite field, i.e. if df.jf has a minimum at M = Mq. In cases 
where d\jf becomes negative, a thermodynamic instabil- 
ity is produced where the magnetization jumps abruptly 
between two values Mi and M 2 at which H m = H m 
[Fig.^a)], determined by a Maxwell construction: H m = 
d M f\ Ml = d M f\ M2 = (f(M 2 )-f{M 1 ))/(M 2 -M 1 ). The 
location of the jump defines a line of first-order transi- 
tions in the H-T plane, which ends at a critical endpoint 
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FIG. 1: (a) Magnetization as a function of the local internal 
magnetic field H m with a jump at field H m . (b) Schematic 
phase diagram in the H-T plane. For n = 0, there is a line of 
first-order transitions terminating at a critical end-point. We 
assume that H m is a decreasing function of temperature, re- 
flecting the fact that entropy increases with increasing magne- 
tization in itinerant metamagnets [f| . The grey area indicates 
the region of phase separation for n > 0. 



[Fig. nTb)]. A quantum critical endpoint occurs if T c is 
suppressed to zero as a function of additional parameters. 

This standard picture must be supplemented to include 
magnetic dipolar interactions. We demonstrate this by 
translating an old argument by Condon |7j into the con- 
text of metamagnetism. The internal magnetic field is 
a combination of the external applied field H cx and the 
field created by the magnetized sample H m « H c * — nM, 
where n is the demagnetization factor depending on sam- 
ple shape and field orientation (n = for a needle-like 
sample oriented parallel to the field and n = 1 for a 
thin film oriented perpendicular to the field). Thermo- 
dynamic stability is lost if H in reaches H m : from the 
low (high) field side, this occurs when H ex — H c i — 
H m + nMi, (H cx = H c2 = H m + nM 2 ). Instead of 
a single first-order transition, there is a finite region, 
H c i < H cx < H c2 where the uniformly magnetized state 
is unstable [Fig. ^b)]. In this region, as we will show, 
magnetization forms domains with M = Mi or M 2 . 

A complete theory is based on the free energy func- 
tional for a spatially varying magnetization M(r) parallel 
to the z-axis |ll| . 



F= d 



f(M) + K (VM)" 



(1) 
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where if is a parameter controlling the stiffness, and 



d 3 r'M(r)A(r-r')M(r'), (2) 



with A(r) = (r 2 — 3z 2 )/r 5 is the dipolar magnetostatic en- 
ergy. The equilibrium magnetization profile M (r), which 
minimizes F — fi f d 3 r H ex M , is determined by a compe- 
tition between domain wall (DW) and dipolar energies. 
Below, we discuss its properties in the different regions 
of the phase diagram. 

Close to the critical temperature, one may expand / 
in powers of m(r) = M(r) — Mo, where Mq is the local 
minimum of d M f. Hence, 

f(M) = /(M ) + H H m + tm 2 + um 4 + 0(m 5 ), (3) 

where fj,oHo — <9a//|j\/ . All the parameters Mo, Hq and t 
are temperature-dependent. Apart from the linear term 
which shifts the critical point to finite fields, we obtain 
the standard theory of an Ising ferromagnet with dipo- 
lar interactions |l2(. There is no third-order term by 
construction, but other odd terms such as m 5 are per- 
mitted. For simplicity, we now consider an infinite film 
of thickness D. For this geometry with the field direc- 
tion perpendicular to the slab, the demagnetization fac- 
tor is n = 1. If the magnetization is independent of 
the z-coordinate, the quadratic part of F in powers of 
m is D Jd 2 qLLj q \mq\ 2 /(2ir) 2 , where m q = Jd 2 r e ,q ' r m(r) 
and Lu q = t + Kq 2 + /i (l - e~ qD )/(2qD), We assume 
that the sample is sufficiently thick that e~ qD may be 
neglected. No matter how small the dipolar interaction 
may be in relation to the spin stiffness, the instability 

1/3 

always occurs at a finite wave vector qo = (2l 2 D) , 
where I 2 = 2K//j,q, giving sinusoidal modulations of the 
magnetization with wavelength d — 2n/qo. On following 
the line H cx = Hq + Mq in the H cx -T phase-diagram 
from high temperatures, there is a second-order transi- 
tion from the uniform to the modulated phase as soon as 

the instability condition t < — |/io (jo) 2 ^ 3 is met. For 
values of the external field different from Hq + Mq, both 
bubble- and stripe-modulated phases may be stabilized 
within mean-field theory and the transitions are weakly 
first-order [J^- In each case, the modulated phase con- 
sists of a smooth, sinusoidal spin-density wave on top of 
the uniform magnetization background. 

When T < T c , the double- well structure of g(M) = 
f(M)—fj,oH m M becomes more pronounced, such that si- 
nusoidal magnetization oscillations are no longer favored. 
The sample splits into domains of constant magnetization 
Mi and M 2 , separated by relatively sharp DWs. The de- 
termination of the domain structures then involves two 
independent problems. The first is the internal structure 
of the DW, which depends on the potential g(M) and the 
stiffness AT, and the second is the global domain pattern, 
which depends solely on the long-ranged dipolar inter- 
action and the characteristic length £ = <r/(/ioAAf 2 ), 



where AM = M^ — M\ is the difference of magnetiza- 
tion of the two phases and o is the DW energy per unit 
surface. The first problem is easily solved. The equilib- 
rium magnetization profile of a single, flat DW satisfies 
dM/dx = y/(g(M) - gJMxffjK, which leads to the fol- 
lowing characteristic DW width and energy: 



A = lim 



f 



M 2 -tAM 



dM 



K 



o ln(l/e) J Ml +eAM y/g{M) - gJMT) 



M, 



2VK I dM y /g(M)-g(M 1 ). 

I Mi 



(4) 



(5) 



The second problem has been studied in great detail in 
the context of uniaxial (Ising) ferro mag nets in the shape 
of thin platelets far below T c 0, If4| . In an interme- 
diate range of the ratio D/£, the phase diagram as a 
function of magnetic field in the sharp-wall limit is very 
similar to that in the critical regime: A stripe phase in 
the center and a bubble phase at the border of the coex- 
istence region separated by first-order transitions. The 
typical domain size is of order d ~ (^D) 1 / 2 , and differs 
from the wavelength 2ir/qQ in the critical regime. The 
domain structure is more complicated for thick samples 
(£ <C D ), where an increasing number of wedge-shaped 
domains are formed close to the sample surface while the 
bulk domain size is typically d - (^D 2 ) 1 / 3 Q. 

In order to interpolate between these two regimes, the 
critical regime close to T c and the sharp-wall regime at 
low temperature, we propose the following Ansatz, which 
depends on three variational parameters: A, R, and d. 
For periodic stripe domains, 



M x ,R,d(x) = Mi + 



7rAAM \ ^ sin(g„ R) cos^x) 

d 4- sinh^ ' l) 



where q v = 2irv/d. For a triangular lattice of cylindrical 
bubble domains, 



M x ,R,d(r) = Mi + 



2ir 2 XRAM ^ Ji(Q V uR) cos(Q v/i • r) 

+ — — — ■ — 



Vzd 2 



sinh 



(7) 

where J\ is a Bessel function of the first kind and Q^^ = 
^-(v, ^^t 1 ^ ) . In both cases d is the period of the struc- 
ture and [i, v are integers. The interpretation of the pa- 
rameters R and A depends on the regime. If A <C R, d, 
the Ansatz consists of domains of radius R with mag- 
netization M2 in a sea of Mi and A corresponds to the 
DW thickness. For A > d, the Ansatz yields sinusoidal 
spin-density waves because higher Fourier modes are sup- 
pressed exponentially. The mean value around which the 
magnetization oscillates is controlled by R and the ampli- 
tude by A. Examples of stripe and bubble configurations 
in both regimes are illustrated in Fig. [21 

After these general considerations, we now consider a 
model for a specific itinerant metamagnet, namely a 2D 
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FIG. 2: (Color online) Mean-field phase diagram as a function of external magnetic field ff ex and temperature T for the model 
of Eq. ||SJ within the Ansatz of @, Q. The parameters chosen were 7 = IOO/10, 5 = — and D = 201, where I = (2Kj fio) 1 ^ 2 ■ 
The insets show the magnetization profile m(r) = M(r) — Mo at four different locations. The transitions are first-order at low 
temperature, weakly first-order in the critical regime T <T C and second-order at T — T c . 



electron system whose Fermi level is close to a logarith- 
mic van Hove singularity (VHS). For this situation, the 
free energy f(M,T) was calculated within Hartree-Fock 
theory in Ref. |||. Here, we present an effective model 
which captures the primary physical features. At the 
Hartree-Fock level, fJ-oX 1 — ^\iS x Pf 1 + Pj 1 ~ 27, 
where pj and pi are the thermally broadened densities of 
states at the Fermi level of up- and down-spin electrons 
and I is the on-site Coulomb energy. At zero temper- 
ature, a singularity of the form d\ [ f ~ l/ln|M — Mq\ 
occurs, where Mq is the magnetization at which one of 
the two electron species has its Fermi level exactly at the 
VHS. The main effect of a small, finite temperature is 
to provide a cutoff to this logarithm. These observations 
motivate the following effective free energy [Eq. QJ]: 



Af(M, T) = (s - 7/ln [(m/m Q f + T/T ] ) 



(8) 



where to and To are positive parameters, m = M — M 
and A/ = / - f\ M=M(] - HaH Q m. Since A/ is an 
even function of m, it follows that H m — Hq for this 
model. The model of Eq. JSJ is meaningful only for 
?7i -C mo and T <C Tq. The parameter 7 is of the or- 
der vW//J,%, where v is the volume of a unit cell and 
W" 1 is the weight of the VHS. The parameter S is pro- 
portional to the doping away from the quantum criti- 
cal endpoint in the mean-field phase diagram neglecting 
Ed 0- Dipolar interactions shift the quantum critical 
endpoint to 5 C — — |/Jo (l/2D) 2 ^ 3 . For 8 > 8 C , there is 
a metamagnetic crossover and for 8 < 8 C a first-order 
metamagnetic transition with a coexistence region below 
T c = T exp {7/(5 — 8 C )} . At zero temperature and as- 
suming sharp DWs A <C d, we find using the Maxwell 
construction and Eqs. Q4I5I8J1 : 



AM 



2m exp{(7 + v / 7 (7-4 ( 5))/(4 ( 5)}, (9) 



X/l = [0.98 + (5/7)] VW/H, (10) 
f/i = [0A2 + O(8/ 1 )](-8)/^y. (11) 

The crossover between the critical and the sharp-wall 
regime occurs at 1 - T/T c w {21/ D) 2 ^ / 8 2 . 

The entire if cx -T phase diagram for the model of 
Eq. (jHJ is shown in Fig. [21 It was obtained numeri- 
cally by using the Ansatz of © and J7J) and minimizing 
F — fio J d 3 r H ex M. Figure |3 shows the mean magnetiza- 
tion as a function of H ex at low temperature. The abrupt 
increase of the magnetization both on entry and on exit 
from the domain phase are clearly visible. Much smaller 
increases appear also at the transitions between bubble 
and stripe phases. The energy difference between the 
bubble and stripe solutions is extremely small, so that the 
precise domain structure will be susceptible to many de- 
tails, e.g. lattice structure, impurities, and sample shape. 
The inset of Fig. [3] shows the temperature dependence 
of the magnetization for constant H cx . Note the increase 
of m with decreasing temperature in the uniform phase 
due to the diverging susceptibility at H m = Hq. Upon 
entering the domain phase, the magnetization increase is 
cut off since the internal field H ln is locked to the nearly 
temperature independent value of H m . A small magneti- 
zation jump occurs at both first-order phase boundaries. 

The existence of magnetic domains is known to in- 
crease the electrical resistivity due to DW scattering. In 
uniaxial ferromagnets, the DW resistivity shows a sharp 
increase upon entering the domain pha se by tuning the 
magnetic field at low temperature [15j. The resistivity 
in the direction perpendicular to oriented DWs is higher 
than in the parallel direction [ljj and even shows a non- 
metallic behavior I n contrast to ferromagnetic do- 
mains, Condon domains differ by the magnitude rather 
than the direction of their magnetization. This implies 



4 



AM(T = 0} 




FIG. 3: Mean magnetization m = M — Mo as a function of the 
external field at T = 0.1 T c . Model and parameters are as in 
Fig. |5| The dashed line indicates the (unstable) solution if the 
magnetization were uniform. The solid line was obtained with 
the Ansatz of Q. Inset: Temperature dependence of the 
magnetization along the line H-H -M = 0.08 AM(T = 0). 



resonance lines as long as the corresponding dynamics is 
slow compared to the experimental time scales. 

In conclusion, we predict that an itinerant electron sys- 
tem in general cannot undergo a first-order phase transi- 
tion without breaking up into magnetic domains within 
a finite range of applied fields. While we cannot rule out 
other mechanisms explaining the unusual feature of the 
low temperature phase in S^R^Oy, our results show 
that a detailed experimental analysis is required to de- 
cide if mechanisms other than Condon domains underlie 
the unusual behavior. 

B.B. would like to thank C. Capan, E. Demler, A.G. 
Green, Z.Q. Mao, B. Normand, C. Pfleiderer, A. Scholl 
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the Swiss National Science Foundation and the Science 
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that properties such as the Fermi velocity vary across 
Condon DWs leading to strong scattering of electrons, 
especially if the Fermi level is near a VHS, where elec- 
tronic properties depend strongly on the band filling. 

Our results are particularly relevant in view of the re- 
cent observation of a new phase in Sr3Ru207 with a field 
oriented along the c-direction The qualitative sim- 

ilarities between the phase diagram (Figs. ^ and [2)) and 
the phase diagram of Ref. Q are striking. Also, if meta- 
magnetism is caused by a logarithmic VHS, then the do- 
main phase is very sensitive to sample purity, because 
the VHS itself is very sensitive to impurities. The au- 
thors of Ref. 0] give a list of five key experimental facts 
which place rather tight constraints on any theory for 
the new phase. We note that all five points are compat- 
ible with magnetic domain formation, (i) DW scatter- 
ing leads to an enhanced resistivity, (ii) the electrons are 
itinerant throughout the phase diagram, (iii) the mean 
magnetization increases both on entry and on exit from 
the domain phase as a function of field (cf. Fig. [31), (iv) 
a rise of the mean magnetization is cut off as the do- 
main phase is entered from high-temperature (cf. inset 
of Fig. I3J), and (v) the domain phase exists only in a 
narrow range of applied fields. From the observed field 
range [i (H c2 - H cl ) w 0.2T and the volume 78A 3 
per Ru, we estimate nAM w 1.3 /is/Ru. While this 
number appears rather large, it is not in error by orders 
of magnitude. Finally we note that we have restricted 
ourselves to the simplest case of a field along the c-axis. 
Other orientations require corrections beyond the scope 
of this letter. /jSR, NMR spectroscopy and micro Hall 
probe technics could allow a direct detection of domain 
formation. For (quasi-) static domains the difference in 
local fields will lead to a characteristic splitting in the 
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